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Abstract 

Local Hamiltonians, the central object of study in condensed matter physics, are the quantum 
analogue of CSPs, and ground states of Hamiltonians are the quantum analogue of satisfying assign- 
ments. The major difference between the two is the existence of multi-particle entanglement in the 
ground state, which introduces a whole new level of difficulty in tackling questions such as quantum 
PCP, quantum analogues of amplification, etc. 

The Lieb-Robinson bound is a sophisticated analytic tool used in condensed matter physics for 
handling quantum correlations in ground states, by bounding the velocity at which disturbances 
propagate in quantum local systems. In this paper we show that the detectability lemma (introduced 
in a different context in Ref. [I]), when viewed from the right perspective, can be used in place of 
the Lieb-Robinson bound for the rich case of frustration free Hamiltonians. The advantage of this is 
that the resulting proofs are simpler and more combinatorial, and may be generalizable to solve some 
of the most fundamental questions in Hamiltonian complexity. Additionally, we give an alternative 
proof of the detectability lemma, which is not only simple and intuitive, but also removes a key 
restriction in the original statement, making it more suitable for this new context. 

Specifically, we use the detectability lemma to give a simpler proof of Hastings' seminal ID area 
law [2] for frustration-free systems. Proving the area law for two and higher dimensions is one 
of the most important open questions in Hamiltonian complexity, and the combinatorial nature of 
the detectability lemma based proof and the resulting simplification holds out hope for a possible 
generalization. We also provide a one page proof of Hastings' proof that the correlations in the 
ground states of gapped Hamiltonians decay exponentially with the distance (once again, restricted 
to frustration-free systems). We argue that the detectability lemma in this form constitutes a basic 
tool for the study of local Hamiltonians and their ground states from a computational point of view. 
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1 Introduction 



Local Hamiltonians and ground states, the central object of study of condensed matter physics, are the 
quantum analogues of the central objects of study in computational complexity: constraint satisfaction 
problems (CSP) and their satisfying assignments. This connection, which ties together two seemingly very 
different areas, is the starting point for the emergence of the new field, Quantum Hamiltonian Complex- 
ity, in which properties of local Hamiltonians and ground states are being studied from a computational 
complexity point of view. Over the past few years, this direction has shed exciting new insights into 
quantum information theory as well as into quantum physics. Of crucial importance here is the differ- 
ence between the quantum and classical domains: the quantum analogue of the satisfying assignment, 
namely the ground state, can exhibit extremely intricate multi-particle entanglement. This additional 
player in the game makes borrowing results from the classical domain to the quantum domain extremely 
challenging, cf. the wide open major open problem of whether a quantum analogue of PCP holds [T|; 
it also opens up completely new directions of research regarding the entanglement properties of ground 
states of local Hamiltonians. 

General quantum states require 2™ complex numbers to describe. One of the major goals of quantum 
Hamiltonian complexity is to derive bounds on the entanglement exhibited in ground states of interesting 
classes of local Hamiltonians; the purpose of those bounds and restrictions on the entanglement is to lead 
to an efficient description and analysis of ground states in cases of interest. There is a beautiful sequence 
of papers using structures called tensor networks, with special cases such as MPS [3J |H [SJ |5] , PEPS [7], 
TN [Sj, and MERA jS], which provide such efficient descriptions in certain cases. 

Area laws constitute one of the most important tools for bounding entanglement in such systems. 
Consider the interaction graph (hypergraph) associated with a local Hamiltonian - it has a vertex for 
each particle and an edge for each term of the Hamiltonian. Intuitively and very roughly, an area law says 
that entanglement is local in this interaction graph in the following sense: consider a subset of particles 
L. Then the entanglement between L and L in the ground state is locally "concentrated" along the edges 
between L and L; more precisely, the area law states that the entanglement entropy across the cut is 
big-Oh of the number of edges crossing between L and L. This is clearly a very strong restriction on the 
entropy, which in the general case would be of order of the number of particles (nodes) in L. Proving 
area laws for typical classes of Hamiltonians is thus a holy grail in quantum Hamiltonian complexity. 

A few years ago, in a seminal paper [2J, Hastings proved that the area law holds for ID systems (i.e., 
when the interaction graph is a path), for gapped Hamiltonian - that is, Hamiltonians whose overall 
spectral gap is of order 0(1). In this case, the area law says that ground state entanglement across any 
contiguous cut is bounded by a constant. From this, one can deduce that the ground state of such systems 
can be described efficiently (by an MPS of polynomial bond dimension - see Ref. [2]). The question of 
whether ground states in two and higher dimensions obey an area law is still wide open. 

Hastings' proof of the ID area law, and many other proofs related to entanglement and correlations 
in ground states, use sophisticated analytic methods. Perhaps the most important of those is the famous 
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Lieb- Robinson bound (LR bound) jTQl QT] , which bounds the velocity at which disturbances propagate in 
quantum local systems; Fourier analysis, and other techniques are important players too. These analytic 
tools constitute a major barrier for a fuller participation by computer scientists in this important aspect 
of Hamiltonian complexity. Also, these analytic techniques seem to inherently involve the dynamics of 
the system in time, according to the Hamiltonian. However, purely from an aesthetics point of view, it 
should be possible to explain kinematic results about the ground state without resorting to dynamical 
arguments (which is what the LR bound is). Or, in other words, without adding the extra dimension 
of time to the problem. In addition, the kinematic problem seems, at least on the surface, to be of a 
combinatorial nature, thereby suggesting a combinatorial solution. 

In this paper we introduce a combinatorial tool to tackle the above mentioned problems, and, in 
particular, to get a handle on correlations and entanglement in ground states of local Hamiltonians. This 
is a simple, basic version of the detectability lemma of Ref. [Tj. We demonstrate that when the system 
is frustration-free, many of the results that rely on the traditional analytic tools can be obtained in a 
much simpler, direct and intuitive way using this tool; we argue that the detectability lemma in this form 
constitutes a basic tool for the study of local Hamiltonians and their ground states from a computational 
point of view. 

Our starting point is the Detectability Lemma (DL) introduced in Ref. [T|. There, the motivation for 
the DL was quite specific: to help translate classical results about CSPs to quantum results about local 
Hamiltonians. It was used to prove a quantum analog of gap amplification (a component of Dinur's proof 
of the PCP theorem [II])- The DL made it possible to sensibly make a statement of the form "If the 
ground state energy is at least k then the probability that it violates at least ck terms of the Hamiltonian is 
bounded below by a constant" . The DL of Ref. [T] holds under the mild assumption (which is essentially 
true in most interesting cases) that each particle participates in a bounded number of terms of the 
Hamiltonian, and therefore the terms of the Hamiltonian can be partitioned into a constant number of 
layers, each consisting of terms acting on disjoint sets of particles. Ref. [T] also required an additional 
technical assumption, that the number of distinct types of terms of the Hamiltonian are bounded. 

Here, we reformulate the DL and put it in a much broader and basic context. Our reformulation of 
the DL asks the following question: consider a gapped frustration-free local Hamiltonian H = Y^tLi Hi 
with < Hi < 1. i.e., the ground energy of H is 0, and the spectral gap is e = 0(1). The frustration- 
free assumption means that the ground state minimizes the energy of every local term, so no term is 
"frustrated". Can we approximate the projection H gs on the ground state, by a "local" operator? 
Such a local approximation would be extremely useful, as it would enable deducing local properties of the 
ground states such as area laws and decay of correlations. Indeed, such an approximation of a projection 
on the ground state is essentially what is done by the traditional analytic tools that use the LR bound, 
as we explain in Sec. 2] The approximation offered by the DL, however, has more of a combinatorial 
flavor, and is therefore much easier to handle. 

A natural first guess of such a local approximation of n ss is the positive semi-definite operator 
G = (1 - ±H), where m is the number of terms in the Hamiltonian. G fixes the ground state, and 



3 



Layer 1 (even) 
Layer 2 (odd) 



Figure 1: Am illustration of a ID system of two- local, nearest neighbors, interactions. The local terms 
(Hi,H2, ■ ■ ■) can be arranged in two layers (even and odd), such that the terms in each layer do not 
overlap. 



shrinks all the orthogonal space to it by a factor; however the shrinkage is very limited, by a factor 
of (1 — e/m). To get a good approximation, one would need to apply this operator polynomially many 
times, and by this we would lose the locality of the operator. Indeed, the expression (1 — H/m) m contains 
products of m overlapping terms whose overall support is of the order the size of the system. Our challenge 
is therefore to get a local operator that preserves the ground state but shrinks the orthogonal subspace 
by a constant factor, rather than by e/m. For simplicity of presentation in the introduction, let us 
consider the simplest scenario, in which the particles are set on a ID chain, and the interactions are two 
local. Denote by Pi the projection on the ground state of the terms Hi. Notice that the terms in the 
Hamiltonian can be partitioned into two layers, the even and odd terms, each acting on disjoint sets of 
particles (see Fig. [T|); 

Denote by H dd the product of the projections on the ground spaces of all odd terms P\, P3, ... and 
TLeven the product 
want. The DL states: 



def 

by H eV en the product for the even terms. Then the operator A = Il ddJieven is the "local" operator we 



Lemma 1.1 ( Detectability Lemma (DL) in ID) Let A d = TL ddJl even , and let %' be the orthogonal 
complement of the ground space. Then 

" A ' M '" 5 WITW ■ (1) 

The DL says that the application of A to any vector moves the vector closer to the ground state of 
H by cutting down the mass in the orthogonal subspace by a constant factor. This implies that n gs , the 
projection into the ground space of H, can be approximated to within exponentially good precision by 
applying the operator A I times: U gs = A e + e~°^\ 

Let us explain why this operator is indeed "local" . When A is applied £ times to some local perturba- 
tion B that acts on the ground state there is a pyramid-shaped "causality cone" of projections that 
is defined by B. These are simply all terms which are graph-connected to the operator B (see Fig. [2|. 
All the projections outside that cone commute with B and can therefore be absorbed in the ground 
state (since Pi\Vi) = jfi)), leaving us with a local operator of support of size 0(£). Effectively, A e acts 
non-trivially only on a region of width 0{t), when applied to B\Q). 

We give here a new simple proof of this reformulation of the DL, in the process dropping the assumption 
of Ref. [I] about the number of distinct types of terms of the Hamiltonian. 



4 




Figure 2: An illustration of the expression A B\fl) in a ID system. B is a local perturbation, applied 
to the ground state |f2). The local terms underneath it correspond to the Pi projections in A e . The 
pink terms are the projections inside the causality cone of B. These terms are graph connected to B, 
and generally do not commute with it. The blue terms are outside the causality cone, and can therefore 
commute with B and be absorbed by 

The proof hinges on the following observation, which we refer to as the norm-energy trade off. Assume 
by contradiction that A does not move a vector \ip), which is orthogonal to the ground state, very much. 
Then A\ip) must be very close to the range of each P i; but since the range of Pi is the null space of the 
local term Hi, this means that the energy (ip\A' HiA\tp) must be small. However, on the other hand, the 
sum of those energies must be larger than e, since A\i/}) is orthogonal to the ground space; this implies 
that the shrinkage must be quite significant, providing an upper bound on the norm of the vector A\ip). 

However, the above argument is not sufficiently strong. Since there are m terms Hi, the energy 
contribution of each term can be as small as e/m; this will lead to a factor of e/m in the lemma, which 
is not strong enough. The key point is that the energy-norm trade off can be applied locally, using the 
tensorial structure of A; we break the movement of to A\tjj) into disjoint sequential steps and then 
relate the contributions to the energy of each of the terms Hi with the shrinkage resulting from each step; 
one might suspect that entanglement could prevent such an analysis in which shrinkage accumulates but 
the point is exactly that the local structure of the problem allows this accumulation to happen. Think 
very simplistically of the state a|00000) + b\ 11111) subjected to the local terms Hi = ll^l^. A projection 
of this state on the ground state of Hi for any one of qubits i, results in a shrinkage by a factor of |a| 2 , 
but once one projection is applied in one location, the shrinkage is exhausted and no more shrinkage is 
to be gained by a projection in another location. That this entanglement related phenomenon does not 
happen in the DL scenario is due to the locality of the operators involved; it highlights that the way the 
state A\tjj) can be entangled is severely limited. 

We demonstrate the applicability of this reformulation of the DL by providing significant simplifica- 
tions of the proof of Hastings' area law in ID [2], using the DL in two key points, bypassing completely 
the analytic methods. By this we hope to make this important result accessible to a wider audience, as 
well as possibly extendable to higher dimensions. The outline of the proof still follows that of Hastings, 
but now becomes much easier to understand; we defer the explanation of how the proof goes and how 
the detectability lemma enters the picture to Sec. [5] 

To give another example, we provide a one page, very simple proof of Hastings' celebrated result that 
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the correlations in the ground states of gapped Hamiltonians decay exponentially with the distance [TT| . 
Unlike the area law, this applies to d-dimensional grids for any constant d. More precisely, consider two 
observables A and B that are local and act on sets of particles that are of distance t on the grid; the 
decay of correlations means that the expectation value of their product is almost as that of the product 
of their expectation, up to an error which decays exponentially in I. 

We mention that at first sight, one might connect the exponential decay of correlations to an intuition 
that entanglement between a region L and its surrounding is "located" only close to the boundary of 
L, and thus scales like the area rather than like the volume. Though an appealing intuition, such an 
implication of exponential decay of correlation to area laws is not known, and indeed quantum expanders 
provide a counter-example to such a naive connection |13| . 

In both of those proofs, the DL replaces a combination of the Lieb- Robinson bound with other analytic 
tools; this works of course only when the DL is applicable, namely, for the rich case of frustration-free 
Hamiltonians. The restriction to frustration-free Hamiltonians may seem quite strong. We note, however, 
that there are various frustration-free systems that are interesting from a physics and a computational 
points of view, such as the ferromagnetic XXZ model, the AKLT model [14], and stabilizer codes such as 
the Toric code |15| . In addition, many of the quantum phenomenon in quantum Hamiltonian complexity 
are revealed already in the context of frustration-free Hamiltonians, and the major open problems in this 
area (e.g., quantum PCP and 2D area law) are wide open already for this case. Much is to be learned 
from studying frustration-free Hamiltonians, before we proceed to the more general case; it seems that the 
simpler combinatorial nature of the DL in this case might provide a new handle to those questions, and 
there are reasons to believe that a proof of an area law for frustration-free systems might be extendable 
to the general case. 

To illustrate how exactly the is DL related to the analytic methods, we start our more technical 
discussion with a toy application comparing the usage of the LR bound to the alternative route offered 
by the DL, in Sec. |U 

Related work and further directions: 

The DL seems to be connected to various diverse scientific areas. The connections to the LR bound 
and other analytic tools used in condensed matter physics are discussed extensively in Sec. |4] one other 
connection is to view of the DL operator A as a special instance of the general Method of Alternating 
Projections (MAP), that was first studied by von Neumann |16j . In that method one applies a fixed 
sequence of projections in order to approach the intersection subspace. In the general setting, the pro- 
jections are not assumed to be local, nor the Hilbert space is assumed to be of finite dimension. In recent 
results j!7) , the convergence rate is given as a function of the Fridriechs angle, which is not easily related 
to a physical quantity. The DL, on the other hand, is a MAP under the special assumption that the 
projections are local, associated with a frustration-free fe-local Hamiltonian, with a convergence rate that 
is given as a function of the spectral gap. It would be interesting to see if more insight can be derived 
from these connections. 



6 



Recently, much attention was given to a quantum algorithm which, given a local Hamiltonian, uses a 
process involving random measurements of the energies of the local terms to approach the ground state 
efficiently (for certain cases) [HI [H] • The algorithm discussed in those papers carries similarities to the 
situation we are handling here, despite the fact that measurements are applied rather than projections, 
and also that the terms are chosen randomly, rather than in some fixed order. It seems that the DL lemma, 
and the energy-norm trade off, could potentially be useful also for the analysis of such algorithms. In 
particular, it would be very interesting to see a version of the detectability lemma which applies for the 
case in which the terms are chosen randomly. 

As discussed above, it is a wide open question to apply the combinatorial tools presented in this paper 
to the major open problems of quantum PCP and area laws in dimensions higher than 1, as well as to 
many other basic open questions in quantum Hamiltonian complexity. 

Paper organization: 

We start with notations and preliminaries in Sec. [U and then proceed to the statement and proof of the 
DL in Sec. [3] In Sec. [4] we provide the example comparing the LR bound approach to the DL one. We 
then proceed to the area law proof in Sec. [SJ and conclude with the one page proof of the exponential 
decay in Sec. [51 

2 Notations and Preliminaries 

We consider a fc-local Hamiltonian H acting on T~L = ((D )®™, the space of n particles of dimension d. 
H = Hi where each Hi is a non-negative and bounded operator that acts non-trivially on a constant 
number of k qubits (hence the term local Hamiltonian). We assume that H has a ground space of 
energy 0, which must therefore also be a common zero eigenspace of all terms Hi. This means that H 
is frustration free. We also assume that H is "gapped", meaning that its lowest eigenvalue is (the 
ground energy) and all the next are equal or larger than some constant e > 0. We denote by H 1 C T~L the 
orthogonal complement ground space of H. Thus T-L' is an invariant subspace for H, and 

H\w > el ■ (2) 

Most of these assumptions, except for perhaps the frustration-free assumption, are very often used in 
condensed matter physics. 

Throughout this paper we further assume that the H^s are projections, and hence would be denoted 
by Qi- We define Pi to be the projection on the ground space of Qi, Pi d = 1 — Qi. The assumption that 
H is made of projections is not actually a restriction because we can reduce any frustration-free, bounded 
and gapped system into that case. Specifically, for H = J2i Hi with < K, and a spectral gap r > 0, 
we first add an appropriate constant to each Hi such that their ground energy is 0. Then for every i we 

def 

define Qi as the projection into the space where the energy of Hi is greater than and Pi = 1 — Qi as 
the projection to the ground space of Hi. Finally, we define the auxiliary Hamiltonian H' = Qi- This 
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system is frustration free because the original ground states would also be ground states in H' with a 
vanishing energy. Moreover, for any state \ip±) £ H' and every Hi, 

and therefore the gap in H' is e > t/K. It follows that all of our results can be applied to bounded 
frustration-free Hamiltonians by replacing the gap e in DL with the scaled version t/K. 

Given a state |</>) and a partition of the qubits to two non intersecting sets, R and L, with corresponding 
Hilbert spaces Hl^Hr, we can consider the Schmidt decomposition of the state along this cut: \<f>) = 
^2jOLj\Lj) ® \Rj)- Here ot\ ~> oti > ... are the Schmidt coefficients. Their squares are equal to the 
non-zero eigenvalues of the reduced density matrices to either side of the cut pL{<t>) and pr(4>), which we 
denote by Ai > A2 > • • • . The Schmidt rank of \<p) is then the number of non-zero eigenvalues A^ (or 
Schmidt coefficients ctj), and the entanglement entropy is the entropy of the set {A,}, or, equivalently, 
the von Neumann entropy of the matrix (</>). A straightforward corollary of the Eckart- Young theorem 
|20j is then that the truncated Schmidt decomposition provides the best approximation to a vector in the 
following sense: 

Fact 2.1 Let \4>) be a vector on ®T-Lr, and let \\ > A2 > . . . be the eigenvalues of its reduced density 
matrix. The largest inner product between \<p) and a norm one vector with Schmidt rank r is Xj . 

3 The detectability lemma: A new proof 

For clarity of presentation, we will prove the DL in the case stated in the introduction: where the particles 
are set on a line and the local terms are two-local involving nearest neighbors. This proof contains all 
the necessary ingredients for the proof of the more general DL in the case where the Hamiltonian has k 
local terms that can be partitioned into g layers; we make the precise statement of the more general case 
at the end of this section. 

We begin with a simple lemma that quantifies the norm-energy trade-off in the simple case of two 
projections X, Y: we show that if the application of XY does not move a vector very much then the 
energy of that vector with respect to 1 — Y must be small: 

Lemma 3.1 Given arbitrary projections X,Y and \v) of norm 1, if\\XYv\\ =l-e then 

\\(l-Y)XYv\\ 2 <e(l~e) . (3) 

The proof is given in the Appendix. Let us now proceed to prove the detectability lemma. 
Proof of Lemma ll.lt 
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Figure 3: 

Suppose £ H' is a norm 1 state that is orthogonal to the ground space, and define \<j>) = A\ip). 
Notice that for every ground state (fi|i4|^>) = and so \<j>) is orthogonal to the ground space. We 
would like to show that 

fi wrw ■ (4) 

We will find both a lower and upper bounds for the energy of \(f>) , (<fi\H\(j>) , which will give us an inequality 
for ||(/>|| 2 , from which Eq. dU will follow. 

The lower bound is straightforward since \<j>) is orthogonal to the ground state, and so 

\<l>) > eUtf . (5) 

We shall now upper bound the energy (<j>\H\<j>) by carefully upper bounding the contributions of 
the individual terms (<fi\Qi\(fi). We begin by noting that these terms are equal to for i odd since 
A = HoddJieven and QiH dd = for any odd i (recall that H even ,H dd are products of the projections 
Pi = 1 — Qi). We now want to bound the contributions coming from the even terms. 

For this purpose we present A in a convenient form, by reordering its terms. We call the triplet 
product of projections (P1P3P2), (P5P7P6), . . . pyramids, and denote them by A, = P^i^P^i-iPa-2] 
The remaining terms are combined to the operator R = P4P& .... See Fig. [3] for an illustration of 
this structure in ID. Notice that by just using the fact that P$ and P, commute when i and j are not 
consecutive, we can write: 

A = A 1 A 2 ...A m R , 

where m is the number of pyramids which is approximately n/4. 

We will use this reordering to bound the energy contribution of the terms Q2, Qe, ■ ■ ■ ; a symmetric 
argument will bound the remaining even terms Q4, Qsi ■ ■ ■ etc. The energy contribution of Qn-2 will be 
related to the amount of movement produced by the A, portion of of the operator .4. 

The key point in providing this bound is this. We view the transformation of \ip) — > \Aij)) = \<p) as a 
series of steps given by the application of the pyramids Aj. Specifically, letting \vi) = AjA^+i • • • A m R\^), 
we consider the transformation \ip) — > R\ip) — > \v m ) — > \v m -i) \vi) = A\ip). The square of the 
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norm of the first state, after applying R, is a, n = Hit^y 2 . Let a, = ||^i|| 2 /||wi+i |j 2 be the "shrinkage" of 
the square norm (or movement) resulting from the application of the ith pyramid, for 1 < i < m. 

We shall now prove, using Lemma T3. 11 that the shrinkage a.i is related to the energy of the operator 
Q at the top of the same pyramid A^: 



, < 1 - Oi ■ 
We write 

MQ4i- 2 |$ = ||(l-P 4 i_2)^|| 2 = ||A 1 .--A. i _ 1 (l-P 4 ,-2)A l «. i+1 || 2 
< ||(l-P 4l _ 2 )A^ +1 || 2 , 

and recall that Aj = P 4 .;_3P 4 i_iP 4 i_2. We can therefore apply Lemma T3. II to (1 — P4j_2)Aj (with 
Y = Pa-2 andX = P 4 i_ 3 P 4 . t _i), and conclude that ||(1 - P 4j _ 2 )A l u i _i || 2 < (l-||Aj?;j_i|| 2 /lki-i|| 2 )IK_i 
(1 — cij)||uj_i|| < (1 — ttj). Consequently (4>\Q4i-2\4 > ) 5: 1 — a i- 

Using this upper bound gives an upper bound for the energy contribution for Qi, i G {2, 6, 10, . . . }: 

with the constraint on the norm \\(j>\\ 2 < O a *- The right hand side is maximized when all the a^s are 

_2_ 

equal to each other, i.e., ai = \\(j)\\ m , and therefore we are left with an upper bound of the energy coming 
from Q 2 + 4- . . . as: 

{4>\{Q2 + Qe + ---)\<t>) <m 



l-U\\ 2/m 



where the last inequality follows from the fac10 that for every x £ [0, 1], we have m [l — .t 1 '™] < ^-yr- 

For the energy of Q 4 + Q$ + . . ., a similar decomposition to A = (P 3 P 5 P4)(P 7 PgP 8 ) ■ ■ ■ (P2P6 •••) can 
be made, thereby upper bounding the energy by 2 1 ^^ . Combining the energy upper and lower bounds 
we therefore get 

e\\6\\ 2 < (</,\H\<t>) < 2 1 ~y , (6) 



and so 

2 

which gives 



'■'*<1-Uf<l-Uf, (7) 



1 This inequality can be easily verified by noticing that f m (x) = m^/x [l — x l ' m ] + x is equal to 1 for x = 1 and has a 
non- negative derivative for x £ [0, 1] and m > 1. 
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The above proof can be easily generalized to other geometries. In the general case, in accordance with 
Sec. [2] we assume we have a fc-local, frustration-free Hamiltonian H = Qi that is made of projections 
and has a spectral gap e > 0. We further assume that each particle participates in a constant number 
of projections, and therefore the Qi can be partitioned into a constant number of g layers; each layer is 
made of projections that do not intersect each other and are therefore commuting. 

Then for each layer we define the projection IT as the product of all P } ■ = 1 — Qj that are in the 
layer, and define the DL operator A by 

A d ^U g ---U 1 . (9) 

Finally, we define f(k,g) to be the number of sets of pyramids that are necessary to estimate the energy 
contribution of all the Qi terms. In the ID case that we proved, we had f(k,g) = 2, because only the 
even layer contributed energy and we needed two sets of pyramids to cover that layer. In the general case 
it is easy to see that f(g, k) can be crudely bounded by f(g, k) < (g — l)k 9 . 
Using the above definitions, the general DL is 

Lemma 3.2 (The detectability lemma) Consider the local Hamiltonian system that is described above. 
Then 

\\A\ W \\<- -Ira- (10) 

[ e //(fc,. 9 ) + l] 1/3 

4 Comparing the Lieb-Robinson bound approach and the de- 
tectability lemma approach 

In this section, we compare the DL with a standard method used in many of the seminal results in 
quantum Hamiltonian complexity, such as Hastings' areas law for ID gapped systems [2], and Hastings' 
exponential decay of correlations proof [TT] . The method combines the use of the Lieb-Robinson bound 
(LR bound), with a Fourier analysis and the existence of a gap, to reveal the locality properties of the 
ground state. More specifically, the method uses these tools to approximate expressions that involve the 
projection operator to the ground state, n gs , by local operators. 

To understand how this is done, let us concentrate on a simple example of locality in the ground state, 
and derive it using both the DL and the LR bound. 

We focus on n gs , the projection on the ground space of H. On the surface, this projector seems very 
far from being local in any sense. Nevertheless, in gapped systems it does possess some locality properties 
that are crucial to the analysis of correlations and entanglement in the ground state. To see this, one 
standartly considers an approximation of H gs by another operator: 

P d ^ -L= [ dt e - t2 ^e- mt . (11) 
V^q J 
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Here q is a free parameter to be chosen as appropriate. For an eigenvector \E) of H with eigenvalue E, 
we have P q \E) = e~ qE ^ 2 \E). Consequently, if the system has a constant spectral gap e > 0, P q indeed 
approximates Tl gs well: 

\\P q -K ga \\ <e-^' 2 . 

We now want to argue regarding the local nature of P q in various contexts. Let us illustrate the 
LR bound approach with a simple example: consider the expression H gs B\Q) « P q B\Q), where |f2) is a 
ground state of the system with zero energy, and B is some local perturbation. It is easy to see that 

P„B\tt) = — L= [ dte- t2 / 2,1 e- im Be lHt \n) = —L= [ dt e"' 2 ' /2q B(t)\Q) . (12) 
\Z2wq J V 27r< 7 J 

B(t) = er lHi ~Be lHt is the time evolution of the perturbation B. The key point is now to use the 
famous LR bound, to approximate it by a "local" operator, i.e., an operator which acts only on the 
"neighborhood" of the particles on which B acts. The following is an immediate corollary of the original 
LR bound, which we omit for sake of brevity. The full statement of the LR bound, together with the 
proof of this corollary can be found in Ref. [21] . 

Theorem 4.1 (Lieb-Robinson bound (LR bound) , adapted [21j ) Given a local Hamiltonian H = 
Hi on n particles, there exists a constant velocity v s.t. B{t) can be approximated by an operator de- 
noted Bi(t) whose support is inside a ball of radius I — vt around the support of B, s.t. 

\\B(t) -B e {t)\\ < \\B\\-e-° w . (13) 

Given a length scale t > 0, we may now set q = £ in Eq. (|T2l) and obtain 

n gs B\n) « p t B\si) pa f dt e - t2/2e B(t)\n) « f dt e - t2/2i B e (t)\n) . 

v2tt£ J\t\<t V27r£ J\t\<t 

In the above series of approximations f=s implies an approximation of up to an error of e~°( e \ The 1st 
approximation follows from the assumption of the constant gap and Eq. (|11[) . The 2nd approximation 
is due to the exponential decay of the filter function e~* 2 / 2 ^, and the 3rd is due to the LR bound. We 
therefore get an exponentially (in £) good approximation to H gs B in the expression H gs B\Q) by an 
operator which is ^-local. 

Let us now derive the same result for the frustration-free case using the DL. First, we approximate 
the ground space projection H gs by applying the DL operator A for m times. By Eq. (fT0|). A leaves the 
ground space invariant while shrinking the orthogonal space by a constant factor. Therefore 

n gs = A m + e-° {m) . (14) 

We now write 

n gs B\n) = A m B\n) + \\b\\ ■ e -° [m) , 
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and consider the expression A m B\VL). By assumption, the system is frustration free, and therefore every 
local projection operator Pi that appears in A leaves |f2) invariant: Pj|fi) = |fi). We now consider the 
"causality cone" of projections in A m that are defined by B. These are simply the projections that are 
graph-connected to B when all the projections in A m are arranged in consecutive gm layers (see Fig. [2]). 
The main observation is that all the projections outside this causality cone commute with B, and can 
therefore be absorbed by |f2). We are therefore left only with the projections of the causality cone, whose 
support size I is proportional to m. In other words, just as in the LR bound method, we found an 
exponentially (in £) good approximation to Tl gs B in the expression Tl gs B\fl) by an operator which is 
£-local. 

This kind of reasoning, with appropriate modifications, is used in both the ID area-law and the 
exponential decay of correlations, that are presented in the following sections. 

5 The area law in ID using the detectability lemma 

Throughout this section, we let H = J] Qi be a 2-local frustration-free ID Hamiltonian that is made of 
projections Qi acting on particles of dimension d. Assume that H has a unique ground state |fi) and a 
spectral gap e, and set d d = 1 — (1 + e/2) -1 / 3 in accordance with the shrinking exponent of the DL, in 
Eq. JT]). We notice that in the interesting limit e — > 0, we have 6 ~ e/6, and generally, for for e € [0, 1], 
have e/8 < 5 < e/6. In order to keep the presentation simple, we shall assume throughout this section 
that e < 1, and prove the following version of a one dimensional area law: 

Theorem 5.1 (Area Law for frustration free Hamiltonians in ID) For any contiguous cut along 
the chain, the entanglement entropy of the ground state \Q) across the cut is bounded by a constant which 
depends on the dimensionality of the particles d and on the spectral gap e; specifically, 

S<^-d 4 / s Qnd) 2 . (15) 
o 

The proof relies on two main lemmas. The first shows that for any cut along the line, there is a 
product state |0i) (g> |02) that has a constant inner product with the ground state: 

Lemma 5.2 (Constant overlap with a product state) For every cut, there is a product state |0i)<8> 
|0 2 ) such that ® (j> 2 \n)\ > fi d = dT e {l - 8Y a/A , with £ d = d 4/5 . 

The second lemma shows that if there exists a product state with a constant overlap with the ground 
state |fi), then |f2) has finite entanglement entropy: 

Lemma 5.3 (Constant overlap with a product state implies finite entropy) If for some cut there 
exists a product state \<p\) <E> \4>2) such that \ {<p\ ® 02|^)| ^ I 1 , then the entanglement entropy of |fi) across 
the cut is bounded by 

S< |(ln4r + 2)lnd . (16) 

o [10 
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Theorem 15.11 then follows easily by combining the two lemmas and using the facts that S > \i and 
In i > 6. We prove Lemma \5. 2 1 in Sec. 15.21 and Lemma [5~3l in Sec. 15.11 

5.1 Constant overlap implies finite entropy (proof of Lemma 15. 3j) 

In this section we prove Lemma l5~3l The DL is clearly the right tool for the task, since it provides a "local" 
operator that can be repeatedly applied to the promised product state \4>i) <g> |</>2) without increasing its 
entanglement rank much, while exponentially decreasing its distance from the ground state. 

The only thing that is not entirely clear is how to get a constant bound on the entanglement entropy 
of the ground state, since a straightforward argument would mean applying the operator non-constant 
number of times to get arbitrarily close to the ground state. The key is to observe that after £ applications 
of the DL we get a state with a bounded Schmidt rank that is close to the ground state, and by Fact 12. 1[ 
this gives us a bound on the sum of the largest d 2i Schmidt coefficients of the ground state. With these 
bounds we can find a pessimistic constant upper bound on the entanglement entropy. We can now proceed 
to the more detailed proof. 

Consider then a cut in the line between the particles io and io + 1, and let Qi a be the local term in H 

def 

that involves io, io + 1- Assume that along that cut, the product state |</>o) = |</»i) <E> |</>2) has a constant 
projection /i on the ground state 

\c/>o)=n\n) + \w) , (17) 

where \w) £ W ' , and ||w|| = (1 — /i 2 ) 1 / 2 < 1. We now apply the operator DL operator A I times on \4>o)- 
We obtain 

where \w^) <E W and ||n>W|| < (1 — 5) e . Let \vg) be the normalized version of A^o)- Then 

IH fi >l > r=k ™ = 1 - - 5)2t) ■ (18) 

VM + (1 - <)) 

This means \ve) are exponentially close to the ground state, as a function of I. 

How entangled are those states? We notice that at each application of A, the entanglement rank of 
the state can only increase by a multiplicative factor of d 2 : for every i ^ io, the projection term Pi in A 
works entirely left to the cut or entirely right to the cut, thereby not increasing the Schmidt rank of the 
state. The only projection in A that may increase the rank is Pj , and as it is a 2-local projection that 
works on (i-dimcnsional particles, it can at most increase it by a factor of d 2 ! 2 ] Consequently, the Schmidt 
rank of \ve) is at most d 2i . 

2 Given a vector \v) = J^j = i \L)j IS \Rj) °f Schmidt rank r and i < io, Pi\v) = Y^—i (Pi\Lj)) CE5 \Rj) and thus has 
Schmidt rank bounded by r (the symmetric argument shows the same result for i > io). For i = io, we can decompose 
Rio = £k=l-^fc ® ^fc with Jfj, acting on the io'th particle and acting on the (io + l)'th particle. Consequently 
Pi \v) = 5Zj =1 i(Xk\Lj)) ® (Yk\Rj)) has Schmidt rank no larger than rd 2 . 
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We have therefore obtained a family of states with Schmidt ranks bounded by d 2f , which are 

closer and closer to |0). Then using Fact 12. ll together with Eq. (fTS|) . it follows that the eigenvalues of 
the reduced density matrix of |f2) along the cut, Ai > A2 > • • ., must satisfy 



1/2 



< \{vt\Sl)\ < (5>i 
J =1 



y/f + (1 - sr 

which implies the following scries of inequalities 

For every £>1: ^ Xj < -^(1-5) 



2C 



j>d 2t + l 



(19) 



(20) 



From here, the desired upper bound on the entropy can be deduced by choosing the distribution of 
maximal entropy which still satisfies the inequalities in Eq. (pH]) . The following lemma, whose proof can 
be found in Appendix [Bl gives one such bound: 

Lemma 5.4 Consider a probability distribution {Xj} whose values are ordered in a non-increasing fash- 
ion, Ai > A2 > . . and let D > 2 be an integer and K > 1,0 < 6 < 1 some constants such that 



for every £>1: ^ Xj < K9 e 

j>D'- + l 



(21) 



Then the entropy of {Xj} is upper bounded by 



S<3 



1*1^ + 1 



ln(l/0) 

Substituting D = d 2 ,6 = (1 — 8) 2 and K = [i~ 2 gives 



]nD 



(22) 



S< 6 



In -5 



^[i-(i-sy] + 
21n[l/(l-5)] 



lnd 



But 1 - (1 - (5) 2 > (5 and ln[l/(l - 5)] > 5, and so 



S < 6 



In 



1 



2r) 



3 1 
lnd< -(In— + 2) lnd 
^0 



where the last inequality follows from the assumption that S < e/6 < 1/6. 



5.2 A product state having constant overlap with (proof of Lemma 15.21) 

The obvious candidate for a tensor product state with a constant overlap with |f2) is the mixed state 
Pl <£> Pr, where pr, is the reduced density matrix of |fi) to the left of the cut, and pr is the reduced 
density matrix to the right. 
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Let us assume for contradiction that the overlap between |fi) and pl ® Pr, and in fact with any tensor 
product state along a certain cut, is less than (1 — <5)^ 4 for some sufficiently large constant I. If the 
overlap is small, then there is a measurement that distinguishes |f2) from p^ <8> Pr with probability of at 
least 1 — (1 — <5) £ / 2 ; this is simply the projection on the ground state, LI gs . 

The challenge is to show that there is a local such measurement, i.e., a measurement confined to a 
local window, which distinguishes these two states almost as well. Using the DL we shall now find such 
local measurement that distinguishes with a slightly worse probability 1 — 2(1 — SY' 2 . 

Let us denote by p l L (respectively p R ) the reduced density matrix of |f2) restricted to the £ particles to 
the left (respectively right) of the cut. Also, let p 2e be p = \ ft)(fl\ restricted to the 21 particles, I on each 
side of the cut. We refer to the state p l L ® p l R as the "disentangled" version of the state p 2e . The following 
lemma shows that under the assumption that |fi) has low overlap with every product state \4>i) C3> Ifo) 
(along a given cut), there exists a measurement confined to the window of 21 particles around the cut, 
that with high probability distinguishes p 2£ from p £ L <£) p R . 

Lemma 5.5 (Existence of a distinguishing measurement) Assuming that the overlap of the ground 
state with any product state satisfies \ {cj>\ ® 4>2\Q) | < (1 — S) 1 ' 4 , there is a measurement that distinguishes 
p 2i from p l L (g> p R with probability 1 — 2(1 — S) e ^ 2 . 

The DL ensures that by applying the layers one by one, we converge to the projection on the ground 
state quickly, and it is this projection that is exactly the distinguishing measurement we want to approx- 
imate. We can thus apply A only t/2 times, approximating the projection on the ground space; now, 
following the intuition explained in the introduction (and in the example of Sec. 3]), only the causality 
cone of the cut should be used in this measurement, and the rest of the operators in those layers are 
swallowed by the state being measured; this amounts to a measurement which is restricted to the —I, i 
interval and still distinguishes well enough. The detailed proof can be found in the appendix. 

The fact that such a measurement exists, distinguishing the original state confined to the 21 window 
from its "disentangled" version, with high probability, must somehow indicate that there is a lot of 
entanglement along the cut, whose disentanglement caused this distinguishability. This can be made 
precise using an information-theoretical argument: 

Lemma 5.6 (Distinguishing measurement implies large difference in entropies) // there is a 
measurement that distinguishes p 2e from p e L ® p R with probability of at least 1 — 2(1 — S) e ^ 2 , then 

s( P i) + s(p R )-s(p 2e )> s -e-i. 

The lemma implies that the entropy in S(p L ) + S(p e R ) is significantly larger than S(p 2e ), implying that 
disentangling along the cut has introduced a lot of new entropy. The proof is simple, based on relative 
entropy; essentially, all it uses is the fact that a measurement that distinguishes with high probability 
two states, implies high relative entropy between the results of the measurements. Once again, details 
can be found in the appendix. 
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To finish the proof of Lemma [5^1 we now need to derive a contradiction. Denote by S(2£) the value 
of S(p 2£ ), with 11 being the segment centered around the cut that provides our contradictory assumption 
(namely, that any tensor product state has less than (1 — (5)^/ 4 inner product with the ground state). 
Under these conditions, Lemma 15.51 applies, and hence also the conditions of Lemma 15.61 apply to this 
segment. Applying Lemma [5.61 we conclude that S(2£) < 2S(£) — |i + 1. We now want to recursively 
apply this inequality, for the £ long segments on both sides of the cut, and then for the £/2-long segments 
within those segments, and so on. The problem is that the cuts now move to different locations within 
the 2£ long window, and so our assumption no longer applies for these cuts. However, if the inner product 
state with any tensor product state is small along the original cut, it can be shown to be quite small 
also along near-by cuts, and so all the above arguments can be applied for those cuts too. This can be 
formalized in the following claim, whose easy proof can once again be found in the appendix: 

Claim 5.7 If \(<j>i <fe|^)| < P" f or oil product states across cut (k, k + 1) then \(xi <8> X2|^)| < fJ>d for 
all product states across any cut (fc + j, k + j + 1) with —£ < j < £. 

We therefore assume by contradiction that the inner product along a given cut is smaller than /x, 
such that iid e = (1 — <5)^/ 4 , and so along all the cuts in the 2£ window we have that the inner product 
of the states is at most (1 — <5) £ / 4 , and hence our assumptions apply. We can therefore use the same 
argument recursively. Since 5(1) < ln(d), we get (for £ a power of 2), S(£) < £ln(d) — ^£\og 2 £ + £ < 
£{ln(d) + 1) - p log 2 £. Choosing £ such that § log 2 £ > In d+ 1 makes S{£) < thus giving a contradiction. 

. / c . / c 2(ln d + 1) 

Using the fact that d > 2, this can be achieved by £ = d >° since d 1 > 2 5 . 

6 Exponential decay of correlations 

Consider now a Hamiltonian H = ^ ■ Qi which is /c-local and set on a d-dimensional grid. Once again, 
we assume that Qi are projections, and that H is frustration free with a unique ground state |fi) (i.e., 
Qi|f2) =0), and a spectral gap e > 0. We wish to show 

Theorem 6.1 Decay of Correlations in ground states of gapped Hamiltonians on a d-Dim 
grid: 

Consider a setting as described above. Let X, Y be two local observables whose distance on the grid 
from each other is m. Denote X = (fl\X\n), Y = (fi|Y|fi). Then 

\(Q\(X - X)(Y - Y)\Q)\ = \(Q\XY\Q) - XY\ < \\X\\ ■ \\Y\\ ■ e~° {m) . (23) 

Proof: Let us now consider two operators: Pj„, P ou t' Pin is defined by applying the DL £ times to Y and 
discarding all projections outside the causality cone of Y. £ is chosen such that the resulting cone will 
not overlap with X (see Fig. [4|. Therefore £ oc to, with the proportionality constant that is a geometrical 
factor. P ou t is the complement of -P;„, i.e., it the layers that one get by applying the DL m times, but 
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Figure 4: An illustration of the statement (f2|Xy|f2) = (il\X P out Pi n Y\£l) . The operator P out is drawn 
in blue color and Pi n is in red. Note that the number of projection layers is proportional to the distance 
between X and Y. 

with a "hole" where the causality cone of Y is. Together, we have P in ■ P out = A 1 - See Fig. 0] for an 
illustration in ID. 

Pim Pout leave the ground-state invariant. In addition, they commute with X and Y respectively, 
hence 

(n\x = (n\xp m , y\q) = PoutYp) , 

and therefore 

(n\xr\n) = (n\xp m p out Y\n) = (n\xA e Y\n) . 

We now recall that A e is in fact an approximation of the ground state projection IT gs (see Eq. (fT4|) in 
Sec. II, 

/ = n 98 + e- w = n 98 + e-°w , 

and so 

\(q\xy\q) - (n\xn gs Y\n)\ < \\x\\ ■ \\y\\ ■ e -° (m) . 

Assuming that the ground state is unique, H gs = |f2)(f2|, and therefore 

\(n\XY\n) - xY\ < \\x\\ ■ \\y\\ ■ e -° (m) . 

■ 
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A Proof of Norm-Energy trade-off, Lemma 13.11 

Proof of Lemma I3.lt 

Set |it;) = Y\v)/\\Yv\\. Then 

-Y)XYv\\ 2 = \\Yvf ■ ||(1 - Y)Xw\\ 2 . 

By definition, \w) is a normalized vector inside the support of Y and therefore for every vector \ip), we 
have ||(1 — Y)ip\\ < ||(1 — \w) (w\)i/j\\ • Plugging this to the equality above, we find 

||(1 - Y)XYv\\ 2 < \\Yv\\ 2 ■ ||(1 - \w)(w\)Xw\\ 2 = \\Yv\\ 2 ■ (w\X(l - \w)(w\)X\w) 
= \\Yv\\ 2 ■ \\Xw\\ 2 ■ (I - \\Xw\\ 2 ) 
= \\XYv\\ 2 ■ (1 - \\Xwf) 
< ||Ari;|| 2 ■ (1 - \\XYvf ) = (1 - e)e , 

where the last inequality follows from the fact that > ■ 

B Upper bound on the Entropy 

In this section we prove Lemma 15.41 which is required to finish the proof of Lemma l5.3l 
Proof: 

Call the set of weights {Xj} for D e + 1 < j < D e+1 the tth block. Then the constraints in Eq. (|20]l 
imply that for every block £ > 1, 

D i+i 

E *i<K#- (24) 

j=D e + l 
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Obviously, by reshuffling the mass within a block we maintain the constraints. Moreover, it is straight 
forward to see that the entropy contribution of every block is maximized when all the weights in it are 
equal. The maximal distribution is therefore a steps function, which satisfies: 

K9 l K 

in block £, Xj < — rn ? = (6 7 'D) e . (25) 

~~ D e+1 — D e D — 1 

We now define £q to be the first block for which KB < -(1 — 9)9: 

In eif l n eK i K | i 

Hlj9j ° lrr(W ~ ln(l/0) + • ( 6) 

We will bound the maximal entropy by bounding the entropy contribution of blocks up to (and 
including) and blocks from £q onwards. The first is easy, as there are D e ° weights in the low blocks: 

5/<41n-D. (27) 

In the high blocks, Xj < -(1 — 0) < 1/e, so we can use the monotonicity of the function —A In A in 
the (0 : 1/e] range to bound the entropy by 

Su < - K&H-j^ziWI ) 1 ] < Yl Ke t bx{D t+1 ) (28) 

£>l £>£ 

K9 e °\nD / 1 A 

= ^9-\!° + —9) (29) 

- lnC ( £o + T^) ' (30) 

where the second equality follows from standard geometric sums identities, and the last inequality follows 



from the definition of Iq. Next, looking at the lower bound of £q in Eq. ([26]) , it takes standard calculus 
to verify that £0 > 1J \l(iJe) + — T^0' an( ^ so — 2-^ohii), and 

S = Si + S n <U \nD . 

Plugging the upper bound of £q from Eq. ([26]) . we get Eq. (|22l . ■ 



C Proof of Existence of Distinguishing Measurement, Lemma 

EH 

Proof of lemma l5.5t 

Let Q = {Qi : Qi acts only on particles in the 2£ interval}. Let LI be a projection onto the ground 
space of all the operators in Q. We will show that {n, 1 — II} is the desired distinguishing measurement. 
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Clearly Tr(Up 2t ) = 1. We would now like to prove that Tr(n^ <g> p R ) is at most 2(1 - <5) £/2 . 
We start by considering, instead of II, the applications of the DL operator A. 

We can write pl ® p R as a convex combination of rank 1 density matrices of product states of the form 
1 0i ) ® \4>2)- By assumption, the overlap with the ground state is \{4>\ <£> </>2|^}| < (1 — <S) £ / 4 . Therefore, 
since we assume a unique ground state, \4>\ <g> 02 ) = c|fi) + (1 — c 2 ) 1 / 2 !^), with c < (1 — <5) £ / 4 and 
perpendicular to the ground space. Then by the DL, applying A for £/2 times, we get Tr(^4 f / 2 |0i)(0i | ® 
|</ ) 2)(02|) < [(1 — <5) £ / 4 ] 2 + (l — <5) f / 2 = 2(1 — <5) £ / 2 , and this remains true when we take convex combinations: 

Tv{A l / 2 p L ®p R ) <2(1-^ 2 . 

Thus Tr(A £ / 2 p) = 1, and Tr(,4^ 2 p L ® < 2(1 - <S) £ / 2 ; this establishes that applying .4 for 1/2 
times distinguishes between p and ph <® Pr with the desired probability. However, we would like the 
measurement to be confined to a short interval. Since we know that Tr(np 2f ) = Tr (Up) = 1, the proof 
will follow from showing that 

Tr(npi ® p l R ) = Tr(n Pi ® p K ) = Tr(lL4 f / 2 ( Pi ® p R )) < Tv(A e / 2 (p L ® . 

The first equality holds since II only acts on the 2£ particles in p L <S) p l R and the last inequality is trivial; 
it is the middle equality which uses the structure of A 1 1 2 . Indeed, let us write A 1 1 2 — AmA^Ar where 
Am = ■ ■ ■ {Pk-2PkPk+2){Pk-iPk+i){Pk) is a "pyramid" of terms centered at the cut A;, and Al and 

S v ' 

£ groups 

A R arc the terms to the left and right of the pyramid respectively, as in Fig. [5] Then UA t l 2 (pi J <E> 
Pr) = UAmAlAr(pl ® Pr)- But since we applied A for exactly £/2 times, every Pj projection in 
Am is also in the 21 window of II. Therefore AmPi = Am and consequently UAm = II. Similarly, 
AlAr(pl ® Pr) = Pl® Pr, and therefore UA 1 / 2 (pl ® Pr) = U(p^ ® Pr), implying the desired equality. 



D Proof of Information Theoretical bound, Lemma 15.61 

Proof of Lemma 15.01 

Let X and Y be {0, 1} random variables that result from applying the measurement II on p and ph ® Pr 
respectively. Then by the Lindblad-Uhlmann theorem [531 H3] , 

S(p{) + S(p R )-S(p 2e ) = S(p 2t \\p{®p R )>S(X\\Y) = ]T Xl \n^. 

i6{0,l} Vl 

In this case, we have X = 1 with probability 1 and Y is 1 with probability a < 2(1 — 5) e / 2 . Thus using 
straight forward analysis, Yli x i In ^ = ln( ^) > ln(l/2) — | ln(l — 5) > |^ — 1, and the result follows. ■ 
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Figure 5: An illustration of the identity HA e / 2 (pL <£> pr) = H(pl ® pn). The left, right and middle sets 
of projections in Fig. (b) are inside the invariant space of pl,Pr. and II respectively. 

E Proof that close cuts behave similarly: Claim 15.7 

Proof of Claim [57ft 

Assume for contradiction that ® is a product state across the cut (k + j,k + j + 1) with j > 
such that |(xi ® X2|^)| > A^- Schmidt-decompose = 2i=i ® where the cut is between 

the first fe particles and the j particles between k + 1 and k + j. By simple algebra, there exists at least 
one i such that \(Li ® i?j ® X2|^)| > /U which violates the hypothesis. ■ 
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